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Abstract 

The  Training  Range  Environmental  Evaluation  and  Characterization 
System  (TREECS™)  was  developed  for  Army  firing  and  training  ranges. 
TREECS™  contains  varying  levels  of  capability  to  forecast  the  fate  and  risk 
of  munitions  constituents  (MC),  such  as  metals  and  high  explosives  (HE), 
located  within  firing  and  training  ranges  and  transported  to  surface  water 
and  groundwater.  The  overall  purpose  of  TREECS™  is  to  provide  environ¬ 
mental  specialists  with  tools  to  assess  the  potential  for  MC  migration  into 
surface  water  and  groundwater  systems  and  to  assess  military  range 
management  strategies  to  ensure  protection  of  human  health  and  the 
environment.  TREECS™  can  be  used  to  assess  best  management  practices 
(BMPs)  for  military  ranges  to  avoid,  reduce,  and  remediate  MC  concentra¬ 
tions  in  receiving  waters.  New  modules  for  simulating  the  effectiveness  of 
BMPs  have  been  developed  for  future  implementation  into  TREECS™.  This 
report  documents  the  final  mathematical  formulations,  model 
implementation  protocols,  and  input  requirements  of  these  modules. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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Acronyms,  Abbreviations,  and  Symbols 


Acronyms  and  Abbreviations 


AOI 

area  of  interest,  such  as  HE  impact  area 

BMP(s) 

best  management  practice(s) 

CASRN 

chemical  abstracts  service  registry  number 

DG 

down-gradient 

diss 

dissolved  MC  mass 

DU 

depleted  uranium 

EL 

Environmental  Laboratory 

EPED 

Environmental  Processes  and  Effects  Division 

ERDC 

Engineer  Research  and  Development  Center 

EQI 

U.S.  Army’s  Environmental  Quality  and  Installations 
Research  Program 

Auxin 

MC  mass  Aux  entering  an  AOI  export  treatment  system 

HE 

high  explosives 

MC 

munitions  constituents,  such  as  metals  and  HE 

part 

particulate  MC  mass 

RDX 

Research  Department  Explosive,  a  high  explosive 

SAFRs 

small  arms  firing  ranges 

SW 

surface  water 

TREECS™  Training  Range  Environmental  Evaluation  and 
Characterization  System 

trp 

File  extension  for  a  TREECS™  project  application  file 

TSS 

total  suspended  solids  concentration 

UI 

user  interface  for  entering  and  viewing  model  inputs 

UXO 

unexploded  ordnance 

VZ 

vadose  zone 

WFF 

water  Aux  file  used  in  TREECS™ 

WQCMB 

Water  Quality  and  Contaminant  Modeling  Branch 

Mathematical  Symbols 

A 

surface  area  of  AOI,  m2 

Ab 

surface  area  of  a  sedimentation  basin,  m2 

BCR 

bioconcentration  ratio,  mg  MC  in  plant/kg  plant  /  mg  MC  in 
soil /kg  soil  (dimensionless),  where  all  mass  is  dry  mass 

Br 

AOI  burn  rate  for  source  removal,  acres/yr 

C 

dissolved  MC  concentration  in  the  reactor  pore  water,  mg/L 

Ci 

dissolved  MC  concentration  in  water  entering  the  reactor, 
mg/L 
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Cl  dissolved  MC  concentrations  in  water  exiting  the  reactor, 

mg/L 

Ct  MC  total  (dissolved  and  particulate)  concentration  in  water 

within  and  exiting  a  sedimentation  basin,  mg/L 
Ct i  total  (dissolved  and  particulate)  influent  MC  concentration 

for  degradation  reactor  and  sedimentation  basin,  mg/L 
Ctt  total  (particulate,  dissolved,  and  vapor)  non-solid  phase  MC 

concentration  within  the  soil  matrix  on  a  total  volume  basis, 
g/m3 

Dx  dispersion  coefficient  in  the  reactor  flow,  m2/day 

E  AOI  soil  erosion  rate,  m/day  or  m/yr 

Fc  combined  surface  water  MC  flux  from  AOI  excluding  soil 

interflow,  g/day  or  g/yr 

Fdis  solid  phase  MC  dissolution  flux,  g/yr 

FdP  factor  used  to  convert  from  total  concentration  on  a  total 

volume  basis  to  dissolved  concentration  in  the  porous  media 
on  a  total  volume  basis,  dimensionless 
Fd  fraction  of  MC  total  concentration  in  water  that  is  in 

dissolved  form,  dimensionless 
Fe  MC  flux  leaving  AOI  due  to  soil  erosion,  g/yr 

Fes  MC  flux  leaving  AOI  due  to  erosion  of  solid  phase  MC 

particles,  g/yr 

Fi  MC  flux  leaving  AOI  due  to  soil  interflow,  g/yr 

Fi  MC  leaching  flux  from  AOI,  g/yr 

Fp  fraction  of  MC  total  concentration  in  water  that  is  in 

particulate  form,  i.e.,  adsorbed  to  TSS,  dimensionless 
Fprecip  precipitation  flux  due  to  dissolved  pore  water  MC 

concentration  within  AOI  exceeding  the  MC  water  solubility 
limit,  g/yr 

Fr  runoff  MC  flux  from  AOI,  g/yr 

Fsbd  dissolved  MC  flux  exiting  a  sedimentation  basin,  g/day  or 

g/yr 

Fsbp  particulate  MC  flux  exiting  a  sedimentation  basin,  g/day  or 

g/yr 

Ftp  influent  and  effluent  MC  flux  of  a  degradation  reactor  that  is 

in  particulate  form  and  is  untreated,  g/day  or  g/yr 
fB  fraction  of  AOI  soil  that  is  burned  each  year,  l/yr 

fn  fraction  of  the  plant  production  rate  that  is  harvested  during 

phytoextraction;  also  fraction  of  AOI  treated  with 
phytoextraction,  dimensionless 


fp  fraction  of  AOI  area  that  is  treated  with 

phytotransformation,  dimensionless 
f's  fraction  of  AOI  soil  that  is  removed  each  year,  l/yr 

fr  fraction  of  plant  uptake  (bioaccumulation)  of  MC  that  is 

transformed,  dimensionless 

ft  fraction  of  total  flow  and  mass  flux  from  AOI  to  surface 

water  that  is  treated  by  a  degradation  reactor  and/or 
sedimentation  basin;  fraction  of  total  flow  and  mass  flux 
from  AOI  to  vadose  zone  that  is  treated  by  a  degradation 
reactor,  dimensionless 

G  plant  production  rate,  kg  dry  plant  mass/m2-yr 

He  Henry’s  law  constant,  atm-m3/g-mole 

Hr  reactor  dimension  (e.g.,  depth)  perpendicular  to  flow  axis,  m 

i  indicator  for  each  year  in  a  time  series 

j  indicator  for  each  MC  in  a  list  of  MC 

Kd  sorption  distribution  coefficient  for  partitioning  an  MC 

between  soil  particles  and  water,  L/kg 
Kdr  sorption  distribution  coefficient  for  partitioning  an  MC 

between  reactor  media  and  water,  L/kg 
Kdw  sorption  distribution  coefficient  for  partitioning  an  MC 

between  water  and  solids  suspended  in  water,  L/kg 
Kh  dimensionless  Henry’s  constant  for  partitioning  between  air 

and  water 

Lr  length  of  the  reactor  along  the  flow  axis,  m 

L(t)  time-varying  solid  phase  MC  mass  loading,  g/yr 

Ms  solid  phase  MC  mass,  g 

Q  water  flow  rate  through  a  degradation  reactor  or 

sedimentation  basin,  m3/day  or  m3/yr 
Qi  soil  interflow  flow  rate,  m3/day  or  m3/yr 

Qi  leaching  flow  through  vadose  zone,  m3/day  or  m3/yr 

Qr  overland  runoff  flow  rate,  m3/day  or  m3/yr 

Qt  total  water  flow  rate  from  AOI  to  surface  water  (excluding 

soil  interflow,  thus,  equals  runoff  flow)  or  to  vadose  zone 
(equals  leaching  flow  Q ),  m3/day  or  m3/yr 
R  reactor  retardation  factor,  dimensionless 

RE  removal  or  trapping  efficiency  of  a  sedimentation  basin  to 

remove  contaminants  from  the  influent,  percent 
Rns  non-solid  phase  MC  source  or  inventory  first-order  removal 

rate,  yr1 
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Rs  MC  solid  phase  source  or  inventory  first-order  removal  rate, 

yr1 

Ru  the  universal  gas  constant  =  8.206  E-5  atm-ms/g-mol-  °K 

SR  MC  mass  selective  removal  rate,  g/yr 

Sr  soil  removal  rate,  metric  tons/yr 

Ta  ambient  absolute  temperature,  °K 

T1/2  half-life  due  to  phytotransformation,  yr 

TSS  total  suspended  solids  concentration  within  and  exiting  a 

sedimentation  basin,  mg/L 

TSSi  total  suspended  solids  concentration  of  the  influent  to  a 

sedimentation  basin,  mg/L 

TSSir  total  suspended  solids  concentration  of  the  influent  to  a 

reactor,  mg/L 
t  time,  yr  or  day 

U  Darcy  transit  speed  through  the  reactor,  m/day 

Vb  volume  of  water  in  a  sedimentation  basin,  m3 

vs  settling  rate  of  suspended  solids  in  a  sedimentation  basin, 

m/day 

Wr  reactor  dimension  (e.g.,  width)  perpendicular  to  flow  axis,  m 

x  distance  along  the  flow  axis  of  a  degradation  reactor,  m 

Zb  surficial  soil  layer  thickness  of  AOI,  m 

Ar  first-order  degradation  rate  of  MC  within  reactor,  day1 

Ad  reactor  degradation  rate  for  dissolved  MC,  day1 

Ap  degradation  rate  for  MC  adsorbed  to  reactor  material,  day1 

v  pore  water  velocity  in  the  reactor  flow,  m/day 

<f>  AOI  soil  porosity,  fraction 

(f>r  porosity  of  the  reactor  media,  fraction 

pb  soil  dry  bulk  density,  g/cnD  or  kg/L,  g/m3,  and  kg/m3 

pbr  dry  bulk  density  of  the  reactor  media,  kg/L 

Pw  density  of  water,  g/cm3 

6w  soil  volumetric  moisture  content  or  ratio  of  water  volume  to 

total  volume,  fraction;  0W  can’t  be  greater  than  soil  porosity 
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1  Introduction 

Background 

The  Training  Range  Environmental  Evaluation  and  Characterization 
System  (TREECS™)  was  developed  for  Army  firing  and  training  ranges. 
TREECS™  contains  varying  levels  of  capability  to  forecast  the  fate  and  risk 
of  munitions  constituents  (MC),  such  as  metals  and  high  explosives  (HE), 
located  within  firing  and  training  ranges  and  transported  to  surface  water 
and  groundwater.  The  overall  purpose  of  TREECS™  is  to  provide 
environmental  specialists  with  tools  to  assess  the  potential  for  MC 
migration  into  surface  water  and  groundwater  systems  and  to  assess 
military  range  management  strategies  to  ensure  protection  of  human 
health  and  the  environment. 

Firing  and  training  range  Best  Management  Practices  (BMPs)  can  be 
implemented  to  reduce  MC  migration  and  media  concentrations.  BMPs 
can  include  range  use  strategies  as  well  as  remediation  alternatives. 
Military  range  BMPs  for  MC  reduction  generally  fall  under  the  following 
five  broad  categories  (Dortch  et  al.  2013): 

1.  Source  loading  management:  reduce  loadings  to  the  area  of  interest  (AOI). 

2.  Source  removal:  remove  MC  mass  from  the  AOI. 

3.  Source  treatment:  treat  MC  mass  within  the  AOI  to  reduce  export  from  the 
AOI. 

4.  AOI  export  treatment:  treat  AOI  MC  export  (transport)  fluxes  as  they  exit 
the  AOI. 

5.  Down-gradient  receiving  water  treatment:  treat  MC  within  down-gradient 
receiving  waters. 

Examples  of  an  AOI  include  the  impact  area  of  fired  medium  and  large 
caliber  munitions,  firing  points,  impact  berms  for  small  arms  firing  ranges 
(SAFRs),  demolition  sites,  grenade  ranges,  etc. 

Source  loading  management  can  involve  rotating  range  use  over  time  for 
example,  or  simply  varying  the  numbers  and/or  types  of  items  fired  each 
year.  Source  removal,  the  second  BMP  category,  includes  soil  excavation 
or  removal,  controlled  burning  on  the  landscape,  phytoextraction,  and 
selective  removal  of  MC  mass.  The  third  BMP  category,  source  treatment, 
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involves  treatment  technologies  within  the  AOI  that  more  fully  sequester 
MCs  in  AOI  soil,  such  as  soil  amendments  that  reduce  metal  solubility  and 
dissolution  (Larson  et  al.  2007),  or  that  enhance  MC  degradation,  such  as 
adding  hydrated  lime  for  alkaline  hydrolysis  of  high  explosives  (Larson  et 
al.  2008  and  Gent  et  al.  2010).  Although  Phytotransformation  is  a  source 
treatment  BMP,  it  can  also  be  classified  as  source  removal.  The  fourth 
BMP  category  is  referred  to  as  AOI  exit  (or  export)  treatment  because  it 
addresses  MC  removal  after  it  leaves  the  AOI  and  before  reaching  down- 
gradient  receiving  waters.  Three  types  of  export  treatment  devices  are 
addressed:  degradation  reactors;  sedimentation  basins;  and  filter  tubes. 
Down-gradient  receiving  water  treatment  refers  to  MC  removal  within  off¬ 
site  waters.  This  fifth  category  can  include  treatment  within  various  types 
of  surface  waters,  such  as  wetlands,  ponds,  lakes,  and  streams,  as  well  as 
groundwater  remediation. 

Objective 

A  previous  report  by  Dortch  et  al.  (2013)  presented  formulations  and 
software  implementation  and  application  guidance  for  modeling  various 
range  BMPs  within  all  five  of  the  above  BMP  categories.  BMP  modules  were 
developed  for  future  implementation  into  TREECS™  following  publication 
of  that  report.  Adaptations  of  the  implementation  recommendations  were 
required.  The  primary  objective  of  this  report  is  to  document  the  final 
mathematical  formulations  and  software  implementation  protocols  of 
TREECS™  BMP  modules.  Much  of  the  coding  for  the  final  modules  had 
been  written  and  tested  at  the  time  of  this  publication,  but  their  final 
implementation  within  TREECS™  had  not  been  completed. 

Report  layout 

Chapter  2  documents  the  final  formulations  and  implementation  protocols 
of  the  new  modules  in  TREECS™  for  evaluating  BMPs.  Chapter  3 
describes  input  requirements  and  output  for  each  BMP  module.  Chapter  4 
provides  a  brief  summary. 
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2  Formulation  and  Implementation  of  BMP 
Modules 

Background 

As  noted  in  Chapter  l,  the  five  general  categories  for  range  BMPs  are: 

1.  Source  loading  management:  reduce  loadings  to  the  AOI . 

2.  Source  removal:  remove  MC  mass  from  the  AOI. 

3.  Source  treatment:  treat  MC  mass  within  the  AOI  to  reduce  export  from  the 
AOI. 

4.  AOI  export  treatment:  treat  AOI  MC  export  fluxes  as  they  exit  the  AOI. 

5.  Down-gradient  receiving  water  treatment:  treat  MC  export  fluxes  within 
down-gradient  receiving  waters. 

Dortch  et  al.  (2013)  explained  that  TREECS™  could  be  used  in  its  original 
form  (prior  to  adding  any  BMP  modules)  to  address  BMP  categories  1,  3, 
and  5.  Although  not  reiterated  here,  Dortch  et  al.  (2013)  also  discussed  the 
appropriate  application  of  TREECS™  to  address  BMPs  within  each  of 
these  three  categories. 

Phytoremediation  includes  phytoextraction,  phytostabilization,  and 
phytotransformation.  Phytoextraction  is  the  uptake  of  MCs  by  plants 
without  significantly  changing  the  MC  chemical  form  and  harvesting  or 
removal  of  the  plants,  thus  extracting  the  MC  from  the  AOI.  Plants  can 
uptake  metals,  which  do  not  transform.  Thus,  phytoextraction  is  a  form  of 
source  removal.  Phytostabilization  is  the  use  of  plants  to  stabilize  the  soil 
to  reduce  erosion.  TREECS™,  in  its  original  form,  can  address  this  type  of 
BMP  by  using  a  different  type  of  land  cover  and  condition  to  estimate  a 
revised  soil  erosion  rate.  Phytotransformation  is  the  uptake  of  MCs  by 
plants  and  transformation  of  the  MC  to  a  different,  less  harmful  chemical. 
Phytotransformation  technically  falls  into  BMP  category  3  or  source 
treatment,  but  it  can  also  be  considered  a  source  removal,  as  is  the  case 
here. 

This  chapter  describes  the  new  modules  required  by  BMP  categories  two 
and  four.  Some  of  the  BMP  formulations  are  the  same  as  presented  by 
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Dortch  et  al.  (2013),  but  they  are  repeated  here  for  completeness  since 
several  formulations  have  changed. 

Source  removal 

Source  removal  can  involve  five  methods:  removal  of  soil  containing  MC; 
burning  of  the  landscape  to  destroy  MC;  phytoextraction,  where  plants 
uptake  MC  and  are  harvested  or  removed  from  the  AOI;  phytotransforma¬ 
tion,  where  plants  uptake  MC  and  transform  it  to  other  forms;  and  selective 
removal  of  MC,  including  unexploded  ordnance  (UXO)  removal.  Soil 
removal  usually  involves  excavating  and  processing  soils  from  impact  berms 
of  SAFRs  to  extract  metals.  Burning  of  the  landscape  has  been  considered 
for  removing  combustible  MCs,  such  as  HE  and  propellants.  Although 
phytotransformation  is  technically  a  source  treatment,  it  is  included  in  the 
source  removal  BMP  category  to  simplify  module  implementation.  Selective 
removal  of  MCs  includes  UXO  dud  removal  and  removal  of  metal  and 
exposed  HE  chunks.  All  five  of  these  removal  methods  can  be  addressed 
with  modifications  to  the  existing  TREECS™  Tier  2  soil  model  as  explained 
below. 

Since  new  loss  terms  are  involved,  all  of  the  BMP  source  removal  modules 
required  changes  to  the  Tier  2  soil  model  that  is  in  TREECS™.  There  are 
three  new  MC  mass  loss  terms  associated  with  source  removal  BMPs. 
These  loss  terms  include:  a  first-order  loss  rate  of  solid  phase  MC  mass;  a 
first-order  loss  rate  of  non-solid  phase  MC  mass  following  dissolution  by 
water;  and  a  zero-order  loss  rate  of  solid  phase  MC  mass.  There  are  two 
versions  of  the  Tier  2  soil  model,  one  for  constant  average,  annual 
hydrology,  and  one  for  time-varying,  daily  hydrology.  Both  versions  were 
modified  to  handle  source  removal  BMPs. 

The  solid  phase  MC  mass  balance  equation  for  AOI  soil  with  zero-order 
and  first-order  losses  due  to  source  removal  BMPs  is 


dMs 

dt 


Ut)-Fdis-Fes+Fprecip 


RSMS  —  SR 


(1) 


where, 

Ms  =  solid  phase  MC  mass,  g 
t  =  time,  yrs 

L(t)  =  time-varying  solid  phase  MC  mass  loading,  g/yr 
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Fdis  =  MC  dissolution  flux,  g/yr 
Fes  =  erosion  flux  of  solid  phase  MC  particles,  g/yr 
Fprecip  =  precipitation  flux  of  MC  due  to  dissolved  pore  water 

concentration  exceeding  the  water  solubility  limit,  g/yr 
Rs  =  first-order  MC  mass  loss  rate  due  to  source  removal  BMPs, 

i/yr 

SR  =  zero-order  MC  loss  rate  due  to  source  removal  BMPs,  g/yr 

The  last  two  terms  on  the  right-hand  side  of  Equation  l  are  two  of  the 
three  new  loss  terms  of  the  AOI  soil  model.  Otherwise,  the  remainder  of 
Equation  l  is  the  same  as  the  original  version  of  the  Tier  2  soil  model  as 
described  by  Dortch  et  al.  (2011).  It  is  noted  that  Ms  is  the  mass  of  MC 
prior  to  dissolution  into  water  where  it  can  partition  among  water,  soil, 
and  air  within  the  soil  matrix. 

The  mass  balance  for  non-solid  phase  MC  mass  with  a  first-order  mass 
loss  rate  due  to  source  removal  BMPs  is  stated  as 


dCtt 

F<ns 

K 

dt 

~  AZh 

% 


■A  F  )  +  Ap 

a  PP )  g  aP 


Ks  ctt 


_ M  p 

ewz„  dp 

F  . 

precip 


(2) 


AZ, 


where  Ctt  is  the  total  concentration  of  non-solid  phase  MC  on  a  total 
volume  basis  (grams/cubic  meter),  and  Rns  is  the  first-order  non-solid 
phase  MC  removal  rate  (per  year)  due  to  source  removal  BMPs.  Dortch  et 
al.  (2011)  defined  all  other  terms  in  the  above  equations;  however,  they  are 
not  repeated  here  for  brevity.  The  only  change  to  Equation  2  from  the 
original  model  (Dortch  et  al.  2011)  is  the  addition  of  the  last  term  within 
the  bracket  on  the  right-hand  side,  Rns- 

Three  new  parameters  have  been  included  within  the  Tier  2  soil  model,  Rs, 
Rns,  and  SR.  The  model  reads  in  these  three  parameters  from  a  new  input 
file  that  is  generated  from  user  inputs  entered  on  the  BMP  input  screens  of 
the  Tier  2  soil  model  user  interface  (UI).  Each  of  these  three  parameters 
can  vary  over  time  and  is  the  sum  of  all  source  removal  BMPs.  The 
calculation  of  the  contribution  of  each  source  removal  BMP  to  these  three 
terms  is  explained  below. 
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Soil  removal 

Depending  on  the  follow-up  treatment,  soil  removal  can  result  in  first- 
order  losses  for  both  solid  phase  and  non-solid  phase  MC  or  only  first- 
order  loss  of  solid  phase  MC.  For  example,  if  soil  is  permanently  removed 
from  the  AOI,  then  both  types  of  losses  occur.  However,  if  soil  is  processed 
on-site  to  remove  solid  MC  particles,  such  as  lead,  and  returned  to  the 
AOI,  then  only  solid  phase  MC  removal  is  involved. 


The  source  removal  rate  Rs  (per  year)  for  solid  phase  MC  due  to  soil 
removal  is  computed  as  follows: 


Rs  ( soil  removal ) 


[Pb  +6wPw)AZb 


f! 


(3) 


where  Sr  is  the  soil  removal  rate  (metric  tons/year),  pi,  is  the  soil  dry  bulk 
density  (grams/cubic  centimeter),  6W  is  the  soil  volumetric  moisture 
content  (a  fraction  that  is  less  than  the  porosity),  pw  is  the  density  of  water 
(approximately  1.0  grams/cubic  centimeter),  and  the  product  AZb  is  the 
AOI  total  soil  volume  (cubic  meters)  including  water  and  air  spaces.  A 
factor  of  1.0E6,  which  is  the  conversion  from  metric  tons  to  grams,  is 
cancelled  by  the  factor  1.0E6,  which  is  required  to  convert  soil  bulk  density 
and  water  density  in  grams/cubic  centimeters  to  units  of  grams/cubic 
meter.  A  is  the  AOI  surface  area  (square  meters),  and  Zb  is  the  AOI  active 
(contaminated)  soil  layer  thickness  (meters).  The  right-hand  side  of 
Equation  3  is  also  equal  to  the  fraction  of  AOI  soil  that  is  removed  each 
year,  //  (per  year). 


The  constraint  placed  on  the  computed  value  of  Rs  (soil  removal)  is  that  it 
cannot  be  greater  than  1.0.  Values  greater  than  1.0  are  not  allowed  since 
this  would  result  in  more  soil  being  removed  in  a  year  than  is  available 
within  the  AOI  active  soil  layer.  The  only  new  input  required  for  Equation 
3  is  Sr,  which  can  vary  over  time.  The  UI  provides  an  input  table  for  pairs 
of  time  (year)  and  Sr  for  that  year.  As  with  other  time-varying  inputs,  at 
least  two  time-pairs  are  required,  and  values  are  held  constant  until  the 
next  update  time.  The  soil  removal  rate  is  not  MC-specific;  thus,  only  a 
single  entry  is  required  for  each  update  year  and  applies  to  all  MCs. 

If  soil  is  permanently  removed  from  the  AOI,  then  the  non-solid  phase  MC 
removal  rate  Rns  is  set  equal  to  Rs  for  this  BMP  with  the  same  constraint  of 
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not  allowing  values  greater  than  l.o.  The  user  declares  whether  Rns  is  zero 
or  equal  to  Rs,  and  this  declaration  applies  to  all  MCs.  The  parameter  SR  is 
zero  for  this  BMP. 

Landscape  burning 

It  is  possible  to  remove  reactive  MCs,  such  as  explosives  and  propellants, 
through  combustion  by  burning  the  landscape.  Burning  is  MC-specific 
since  some  MCs,  such  as  metals,  do  not  burn.  The  solid  phase  MC  removal 
rate  due  to  burning  is  computed  from 


R  ( burning )  =  4,047  =  f'  (4) 

where  Br  is  the  burn  rate  (acres/year),  and  the  factor  4,047  is  the 
conversion  from  acres  to  square  meters.  The  value  of  Rns  is  set  to  equal  to 
Rs  for  burning.  Computed  values  of  Rs  (as  well  as  Rns )  for  burning  are  not 
allowed  to  be  greater  than  1.0  since  Rs  is  equivalent  to  the  fraction  of  the 
AOI  landscape  that  is  burned  each  year,  fB  (per  year).  The  user  must  select 

which  MCs  are  burned  since  some  MCs,  such  as  metals,  do  not  burn.  An 
input  table  is  provided  for  entering  pairs  of  time  and  burn  rate.  The  value 
of  SR  is  zero  for  this  BMP. 

Phytoextraction 

It  is  assumed  that  plants  can  take  up  only  dissolved  MCs;  therefore,  the 
value  of  Rs  is  set  to  zero  for  this  BMP.  The  value  of  SR  is  also  zero  for  this 
BMP.  The  value  of  Rns  is  computed  as  follows  to  account  for 
phytoextraction: 


Rns  (phytoextraction) 


Fdp  fH  G  BCR 
zb  Pb 


(5) 


where 

Fdp  =  dimensionless  factor  to  account  for  the  soil  pore-water- 
dissolved  portion  of  the  total  non-solid  phase  MC 
concentration  in  soil 

fn  =  fraction  of  the  plant  production  rate  that  is  harvested  each 
year  (also  fraction  of  AOI  used  for  phytoextraction) 

G  =  plant  production  rate,  kg  dry  plant  mass/m2-yr 
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BCR  =  bioconcentration  ratio,  mg  MC  in  plant/kg  plant  /  mg  MC  in 
soil/kg  soil  (dimensionless),  where  all  mass  is  dry  mass 

Other  terms  are  the  same  as  defined  previously.  Assuming  that  all  plants 
are  harvested  each  year  within  an  area  being  treated  with  phytoextraction, 
the  fraction  of  the  plant  production  rate  that  is  harvested,///,  is  the  same  as 
the  fraction  of  the  AOI  that  is  treated  with  phytoextraction;  thus,  this 
fraction  is  not  allowed  to  be  greater  than  l.o. 

The  formulation  for  computing  FdP  is 

O 

7/  _  _ UW _  /0\ 

dP~Sw+(<f-6w)KH+pbK<I 


where 

6w  =  soil  volumetric  moisture  content  or  ratio  of  water  volume  to 
total  volume,  fraction;  6W  can’t  be  greater  than  soil  porosity 
<f)  =  sediment  porosity  or  ratio  of  void  volume  to  total  volume, 
fraction 

Kd  =  sorption  distribution  coefficient  for  partitioning  a  constituent 
between  soil  particles  and  water,  L/kg 

Kh  =  dimensionless  Henry’s  constant  for  partitioning  between  air 
and  water 


All  variables  in  Equation  6  are  known  by  the  UI  for  computing  FdP.  Kh  is 
computed  by  the  soil  model  UI  from  Henry’s  law  constant  He 
(atmosphere-cubic  meter/gram  molecular  weight),  which  is  supplied  by 
either  the  MC  data  base  used  in  TREECS™  or  the  user: 


Kh 


He 

RuTa 


(7) 


where  Ru  is  the  universal  gas  constant  =  8.206  E-5  atm-m3/g-mol-  °K,  and 
Ta  is  ambient  absolute  temperature  (degrees  kelvin),  or  soil-water  matrix 
temperature  in  this  case,  which  is  the  matrix  temperature  in  degrees 
Celsius  plus  273.  All  of  the  right-hand  side  parameters  in  Equations  5-7 
are  presently  soil  model  UI  inputs  or  are  computed  from  inputs  with  the 
exception  of//,  G,  and  BCR,  which  are  all  MC-dependent.  The  first  two 
parameters  indirectly  depend  on  MC  due  to  different  plants  being  used  for 
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different  MC,  and  BCR  is  directly  MC  and  plant  dependent.  It  is  noted  that 
FdP  can  vary  over  time  when  time-varying  hydrology  is  used  due  to  time- 
varying  soil  water  content.  Soil  water  content  is  available  within  a  time- 
varying  hydrology  input  file  used  by  the  soil  model. 

The  fraction  of  the  AOI  used  for  phytoextraction  can  vary  by  MC  and  over 
time.  Thus,  an  input  table  similar  to  before  is  available  to  enter  time 
(years)  and /h  input  pairs  for  each  MC.  A  check  is  made  to  ensure  that 
none  of  thefH  inputs  for  any  MC  is  greater  than  1.0  for  all  time  pairs.  The 
assumption  is  made  that  G  and  BCR  are  constant  over  time  for  each  MC. 
The  values  for  G  and  BCR  are  not  assumed  to  be  the  same  as  those  entered 
for  phytotransformation. 

Burning  for  phytoextraction  is  assumed  to  be  independent  of  burning  the 
landscape  for  source  removal.  The  method  of  harvesting  or  plant 
extraction  is  not  a  model  feature.  Similar  to  burning,  phytoextraction  may 
not  apply  to  all  MCs;  thus,  the  user  is  required  to  select  which  MCs  are 
applicable  to  phytoextraction. 

Phytotransformation 

Phytotransformation  is  a  BMP  that  is  MC-dependent  and  can  vary  over 
time.  Phytotransformation  is  assumed  to  occur  only  for  pore-water, 
dissolved  MC;  thus,  the  solid  phase  MC  removal  rate  Rs  is  zero  for  this 
BMP.  The  value  of  SR  is  also  zero  for  this  BMP.  The  non-solid  MC  removal 
rate  Rns  for  this  BMP  is  computed  from 


Rns  (phytotransformation)  =  fp 


Fdp  G  BCR  fT 
zb  Pb 


(8) 


where 

fp  =  fraction  of  AOI  area  that  is  treated  with  phytotransformation, 
dimensionless 

fr  =  fraction  of  plant  uptake  (bioaccumulation)  of  MC  that  is 
transformed 

All  other  terms  have  been  defined  preciously,  but  it  is  emphasized  that 
values  used  for  G  and  BCR  in  phytotransformation  can  be  different  from 
those  used  in  phytoextraction.  The  soil  UI  already  had  inputs  for  Zb  and  pb, 
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and  information  is  available  for  computing  FdP,  as  explained  above,  for 
phytoextraction;  thus,  the  input  screen  for  phytotransformation  includes 
inputs  for  fp,  G,  BCR,  andfr.  Each  of  these  four  inputs  can  be  MC-specific. 
Plant  production  rate  can  be  independent  of  MC,  but  different  plants  could 
be  used  for  different  MC;  thus,  G  could  depend  indirectly  on  MC. 

Phytotransformation  can  vary  with  time  throughout  the  simulation.  The 
inputs  G,  BCR,  and/r  are  assumed  to  be  constant  over  time  for  each  MC. 
However,  the  input  fp  is  allowed  to  be  time  varying  for  each  MC;  thus,  time 
pairs  for^p  must  be  entered  for  each  MC.  A  check  is  included  in  the  UI  to 
ensure  that  none  of  the  fp  inputs  for  any  MC  is  greater  than  1.0  for  all  time 
pairs  since  treating  more  area  than  the  AOI  is  not  allowed.  Additionally, 
the  fraction/r cannot  exceed  l.o. 


The  phytotransformation  rate  Rns  is  similar  to  a  degradation  rate.  The  half 
life  T1/2  (years)  for  this  degradation  rate  can  be  obtained  from 


T  = 

Jl/2 


0.693 

R 


(9) 


The  time  versus^  table  displays  the  computed  removal  rate,  Rns,  and  half- 
life  for  phytotransformation. 

Selective  MC  removal 

Selective  MC  removal  includes  UXO  duds  as  well  as  chunks  of  exposed 
metals  and  HE.  Duds  are  munitions  that  do  not  encounter  any  explosion. 
Sometimes  range  clearance  personnel  remove  duds  from  firing  ranges. 
Although  UXO  removal  is  performed  primarily  to  remove  intact  duds, 
there  can  be  exposed  chunks  of  HE  from  low  order  detonations  that  can  be 
removed.  The  extraction  of  metal  particles  from  impact  berms  of  SAFRs 
via  soil  processing  can  be  addressed  as  soil  removal,  but  it  can  also  be 
addressed  via  the  Selective  MC  Removal  of  Chunks  BMP.  This  BMP 
module  was  designed  to  allow  selective  removal  of  MC  mass,  and  is 
handled  with  the  zero  order  MC  mass  removal  rate  SR  (grams/year). 
Although  metal  particles  such  as  lead  in  SAFRs  can  be  addressed  with  this 
BMP,  it  was  really  intended  for  exposed  chunks  of  HE  and  large  metal 
chunks  typically  associated  with  spotting  projectiles  that  do  not  explode, 
such  as  those  containing  depleted  uranium  (DU)  that  were  used  in  the 
Davy  Crockett  training  round. 
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Presently,  duds  only  contribute  to  MC  residue  loading  if  percentages  for 
sympathetic  duds  and  their  yield  are  set  greater  than  zero  in  the 
Operational  Inputs  screen  on  the  Site  Conditions  tab.  TREECS™  presently 
does  not  handle  corrosion  and  perforation  of  duds  that  can  expose  HE. 
Sympathetic  duds  are  low  order  detonations  due  to  another  round 
exploding  close  enough  to  the  dud  to  cause  it  to  explode  or  causing  the  dud 
casing  to  open,  thus  exposing  HE.  Although  the  default  values  are  zero,  the 
user  can  enter  the  percentages  for  sympathetic  duds  and  their  yield.  If  the 
percentage  of  sympathetic  duds  remains  zero,  this  can  mean  that  either  all 
of  the  duds  are  removed  or  none  of  them  exploded  sympathetically; 
therefore,  HE  source  removal  associated  with  UXO  dud  clearance  does  not 
require  any  modifications  to  TREECS™.  The  user  can  simply  leave  the 
percentage  of  sympathetic  duds  set  to  zero,  which  has  the  same  effect  as 
removing  all  duds. 

Selective  removal  of  exposed  HE  and  metal  chunks  can  be  addressed  by 
the  user  specifying  which  MCs  are  selectively  removed  and  setting  values 
for  SR,  i.e.,  the  amount  of  solid  phase  MC  mass  that  is  removed  from  the 
AOI  each  year  for  each  MC.  Thus,  SR  is  MC-specific,  and  a  table  is 
provided  for  entering  yearly  time-SR-pairs. 


Implementation 

Each  of  the  five  source  removal  BMPs  are  implemented  as  options  in  the 
Tier  2  soil  model  UI  under  the  Source  Removal  BMPs  tab.  It  is  possible  to 
include  any  and  all  of  these  five  types  (i.e.,  soil  removal,  burning, 
phytoextraction,  phytotransformation,  and  selective  MC  removal  of 
chunks).  However,  the  first  four  source  removal  types  are  assumed  to  be 
mutually  exclusive.  Mutual  exclusion  means  that  if  a  portion  of  the  AOI  is 
treated  with  one  method,  such  as  soil  removal,  then  that  area  cannot  be 
treated  with  any  other  method,  such  as  burning,  phytoextraction,  or 
phytotransformation.  There  is  an  overall  criteria  check  in  the  soil  model 
UI  to  ensure  that  the  user  is  not  using  more  area  for  treatment  than  exists 
within  the  AOI.  The  UI  for  soil  removal  or  burning  adds  the  two  removal 
rates  and  fB  to  ensure  that  their  sum  does  not  exceed  1.0  for  all  time 

pairs.  Furthermore,  the  sum  of  AOI  treatment  fractions  (and/or  fraction 
rates)  for  the  first  four  removal  methods  should  not  exceed  1.0  for  all 
years.  The  overall  criteria  can  be  expressed  mathematically  as  follows, 


M  r  M  r  r'  r'  ^  .  _ 

maX  j— 1  fp, j  +  max  y=i  fntj  +  /Si  +  fBi  <1-0 


(10) 
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where  j  represents  each  MC,  M  is  the  total  number  of  MC,  and  i  represents 
a  specific  year  with  input.  Equation  to  must  be  satisfied  for  each  input 
year,  and  it  conceptually  states  that  if  a  portion  of  the  AOI  is  treated  with 
one  method,  it  cannot  be  treated  with  any  other  method  in  a  given  year. 
The  purpose  of  Equation  to  is  to  avoid  over  counting  for  different 
treatment  options.  Equation  to  is  not  invoked  until  after  all  inputs  have 
been  entered  and  the  user  saves  the  soil  model  UI  inputs.  If  Equation  to  is 
violated,  the  UI  displays  a  screen  stating  “Source  removal  BMPs  should 
not  be  applied  for  an  area  larger  than  the  AOI  in  a  given  year.  Inputs 
should  be  adjusted  to  reduce  BMP  treatments  to  satisfy  this  constraint. 
Prior  inputs  have  not  been  saved.”  At  this  point,  the  user  must  adjust 
inputs  to  satisfy  the  constraints  in  order  to  save  inputs. 

It  is  noted  that  the  implementation  of  BMP  treatment  options  results  in 
conservative  MC  removal  rates  since  the  AOI  is  treated  as  a  single 
homogeneous  area.  For  example,  if  half  of  the  AOI  soil  is  removed  each 
year,  then  all  of  the  MC  source  mass  should  be  removed  after  two  years 
assuming  that  there  is  no  additional  source  of  MC  and  only  the  contami¬ 
nated  half  is  removed  the  second  year.  However,  with  a  homogenous  AOI 
without  spatial  discretization,  61  %  of  the  MC  mass  will  remain  after  the 
first  year,  and  37  %  will  remain  after  two  years.  The  lumped  spatial  aspect  of 
this  model  results  in  a  slower,  exponential  loss  in  MC  mass. 

All  five  source  removal  BMP  options  involve  inputs  as  time  pairs.  All  of  the 
input  pairs  are  MC  dependent  except  for  soil  removal.  The  Tier  2  soil  model 
UI  computes  values  for  Rs  and  Rns  for  each  BMP  option,  for  each  MC,  and 
for  each  update  year.  The  UI  then  sums  each  of  these  two  parameters  over 
all  BMPs  for  each  MC  and  each  update  year  and  writes  the  sums  to  a  file 
that  the  Tier  2  soil  model  subsequently  uses.  Additionally,  SR  values  are 
written  to  this  file  for  each  MC  and  each  update  year.  The  new  file  consists 
of  two  header  lines  for  a  title  and  other  case  information  followed  by  a  line 
with  the  MC  name,  its  chemical  abstract  surface  registry  number  (CASRN), 
and  the  number  of  lines  in  the  time  series.  Each  line  of  the  time  series  has 
the  values  for  the  year,  Rs,  Rns,  and  SR.  Figure  1  shows  an  example  of  the 
new  file. 
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Figure  1.  New  input  file  for  Tier  2  soil  model  with  source  removal  BMPs. 


Example  test  case  for  BMPs 

Data  includes  year,  Rs(l/yr),  Rns(l/yr),  and  SR(g/yr)  for  each  constituent 

lead, 7439921, 2 

0.0,0.1,0.1,2500. 

100.0,0.1,0.1,2500. 

RDX, 121824, 2 
0.0,0.2,0.1,1400. 

100.0,0.2,0.1,1400.| 


AOI  export  treatment 

AOI  export  treatment  involves  the  remediation  of  MC  after  it  exits  (i.e., 
exported  from)  the  AOI  in  transit  to  down-gradient  receiving  waters.  AOI 
export  methods  include  surface  water  runoff  from  the  AOI,  including  soil 
erosion,  and  leaching  from  AOI  soil  to  the  vadose  zone.  Mathematical 
models  were  developed  as  described  by  Dortch  et  al.  (2013)  for  assessing 
two  AOI  export  treatment  options,  a  degradation  reactor  and  a  sedimenta¬ 
tion  basin.  The  degradation  reactor  can  be  used  to  treat  either  surface 
water  export  or  export  to  the  vadose  zone.  Additionally,  the  two  treatment 
methods  can  be  used  in  tandem  (sedimentation  followed  by  reactor)  for 
surface  water.  The  formulations  for  these  two  BMP  options  are  presented 
below.  The  formulation  for  the  reactor  has  not  changed  from  the  one 
presented  by  Dortch  et  al.  (2013),  but  it  is  repeated  here  for  completeness. 
The  formulation  for  the  sedimentation  basin  was  revised  from  the  original 
formulation  by  Dortch  et  al.  (2013). 

There  is  a  third  AOI  export  treatment  option  that  can  be  modeled,  which 
involves  using  filter  tubes  to  trap  or  remove  particulate  and  dissolved  MC 
fluxes  as  they  exit  the  AOI.  Sand  filter  tubes  are  an  effective  means  of 
trapping  particulate  metals  (i.e.,  metals  adsorbed  to  suspended  solids). 
Sand  mixed  with  other  materials  can  remove  particulate  and  dissolved 
MC.  In  most  cases,  filter  tubes  remove  practically  all  exiting  particulate 
MC  and  dissolved  MC  (with  mixed  filter  material).  Thus,  assuming  total 
MC  removal  with  filter  tubes  is  reasonable.  With  this  assumption,  no  new 
TREECS™  modules  are  required  to  include  the  filter  tube  export 
treatment  BMP.  This  BMP  can  be  modeled  by  using  the  User  Defined  WFF 
module  within  the  Advanced  Tier  2  option  of  TREECS™  as  described  by 
Dortch  (2014).  This  module  is  used  in  place  of  the  Tier  2  soil  model  and 
allows  the  user  to  specify  MC  dissolved  and  particulate  fluxes  leaving  the 
AOI.  The  Tier  2  soil  model  outputs  dissolved  and  particulate  fluxes  leaving 
the  AOI  without  filter  tubes.  Users  can  retrieve  those  outputs  from  a 
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previous  run  without  BMPs  and  modify  them  as  inputs  for  the  User 
Defined  WFF  module  as  used  for  the  filter  tube  BMP. 

Degradation  reactor 

This  model  accounts  for  dissolved  MC  removal  along  the  flow  axis  of  the 
reactor.  The  reactor  can  be  oriented  horizontally  or  vertically  to  treat  AOI 
runoff  or  infiltration,  respectively.  The  reactor  is  assumed  to  be  a  water- 
saturated,  porous  medium  with  uniform  dimensions  perpendicular  to  the 
flow  axis  (i.e.,  width,  Wr  (meters),  and  depth,  Hr  (meters),  for  horizontal 
flow)  and  uniform  media  properties,  such  as  reactor  porosity,  </>r.  MC 
concentration  is  assumed  to  be  uniform  except  along  the  flow  axis. 

From  Zheng  and  Wang  (1999),  the  time-varying,  one-dimensional  reactive 
transport  equation  including  sorption  onto  the  solid  media  and  degradation 
is 


„0C  8C 

R - bv - 

dt  dx 


D, 


d2c 

dx2 


-xdc-x„PbrKdr  c 


4 


(ii) 


where 

C  =  dissolved  MC  concentration  in  the  reactor  pore  water  (mg/L) 
t  =  time  (days) 

x  =  distance  along  the  flow  axis  of  the  reactor  (m) 

(fh-  =  porosity  of  the  reactor  media,  void  or  water  volume  /  total 
volume 

Dx  =  dispersion  coefficient  in  the  reactor  flow  (m2/day) 

R  =  reactor  retardation  factor  due  to  sorption  (dimensionless) 
v  =  pore-water  velocity  of  the  reactor  flow,  m/day;  v  =  U/0r,  where 
U  is  the  Darcy  transit  speed  (m/day) 
pbr  =  dry  bulk  density  of  the  reactor  media  (kg/L) 

Kdr  =  distribution  coefficient  for  sorption  partitioning  in  the  reactor 
(L/kg) 

Ad  =  degradation  rate  for  dissolved  MC  (day1) 

Ap  =  degradation  rate  for  MC  adsorbed  to  reactor  material  (day1) 

For  a  given  water  inflow  rate  Q  (cubic  meters/day),  the  Darcy  transit  speed 
U  (meters/day)  through  the  reactor  can  be  computed  from  U  =  Q/(WrHr). 
Equation  11  assumes  equilibrium  partitioning  and  first-order  reaction 
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kinetics.  With  linear  partitioning,  the  reactor  retardation  factor  is 
expressed  as, 


R  =  1 


Pbr^dr 

i 


(12) 


Neglecting  dispersion  and  assuming  that  the  degradation  rates  for 
dissolved  and  adsorbed  MC  are  the  same  within  the  reactor  (Ad  =  Ap  =  Ar), 
Equation  n  can  be  written  as 


dC  v  DC 
dt  R  dx 


(13) 


Assuming  conditions  are  at  steady-state,  Equation  13  reduces  to  the 
following  analytical  solution  for  x  =  Lr,  which  is  the  total  length  of  the 
reactor  along  the  flow  axis, 


CL  =  C,.  exp 


ARL 


(14) 


where  Ci  and  Cl  are  the  dissolved  MC  concentrations  in  water  entering  and 
exiting  the  reactor,  respectively.  This  rather  simple  analytical  solution  can 
be  applied  over  time  with  time-varying  updates  for  flow  rate  and  entering 
MC  concentration  (loading),  but  the  output  concentration  is  assumed  to  be 
at  steady-state  with  respect  to  each  input  update,  which  means  that 
steady-state  is  reached  rather  quickly  relative  to  the  flow/loading  updates. 
This  is  expected  to  be  the  case  for  most  of  the  time,  especially  if  annual 
loadings  are  used.  This  model  does  not  address  the  fate  of  particulate  MC 
entering  the  reactor.  Particulates  could  potentially  settle  in  the  reactor  and 
degrade,  but  this  feature  is  not  supported.  If  substantial  quantities  of 
particulate  MC  are  expected,  then  consideration  should  be  given  for  either 
a  sedimentation  basin  or  tandem  treatment  with  sedimentation  followed 
by  a  reactor. 

Sedimentation  Basin 

The  purpose  of  this  model  is  to  predict  the  trapping  and  resulting  exit 
concentrations  of  MCs  that  are  passed  to  a  sedimentation  basin  via 
overland  flow.  MCs  with  relatively  high-suspended  sediment  -  water 
partitioning  distribution  coefficients  ( Kdw ),  such  as  metals,  would  be  most 
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appropriate  for  this  type  of  treatment.  The  formulation  of  this  model  is 
described  below. 

The  simplifying  model  assumption  that  is  made  is  the  sedimentation  basin 
is  fully  mixed,  i.e.,  a  batch  reactor.  The  only  removal  process  is  settling  of 
particulate  MC  (i.e.,  adsorbed  to  total  suspended  solids,  TSS).  The  mass 
balance  of  MC  within  the  sedimentation  basin  is  stated  as  follows, 

dC 

Vb~jf  =  QCn  -  QCt  -  vsAbFpCT  (15) 


where 

Ct  =  MC  total  (dissolved  and  particulate)  concentration  in  water 
within  and  exiting  the  basin  (mg/L) 

Cn  =  MC  total  concentration  in  water  entering  basin  (mg/L) 
t  =  time  (days) 

Fp  =  fraction  of  MC  total  concentration  in  water  that  is  adsorbed  to 
TSS 

Q  =  water  flow  rate  into  and  exiting  the  basin  (m3/day) 

Ab  =  water  surface  area  of  the  basin,  or  volume/depth  (m2) 

Vb  =  water  volume  of  the  basin  (m3) 

vs  =  settling  rate  of  suspended  solids  (m/day) 

The  fraction  of  MC  total  concentration  that  is  particulate  Fp  is  computed 
from 


10  6  TSS  Kdw 
1  +  10  6  TSS  Kdw 


(16) 


where  TSS  is  the  total  suspended  solids  concentration  in  milligrams/liter, 
and  Kdw  is  the  water-TSS  distribution  coefficient  (liters/kilogram)  for 
sorption  partitioning  of  dissolved  MC  to  TSS.  The  factors  io-6  in  Equation 
16  are  the  conversion  from  mg  to  kg  that  is  required  to  convert  TSS  in 
mg/L  to  kg/L. 

A  TSS  mass  balance  must  be  conducted  for  the  sedimentation  basin  using 
an  analogous  form  of  Equation  15  where  Ct  and  Cn  in  Equation  15  are 
replaced  with  TSS  and  TSSi,  respectively,  and  Fp  is  set  to  1.0.  TSSi  is  the 
basin  influent  TSS  concentration  (mg/L  and  g/m3)  and  is  computed  from 
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TSSi 


1.0  E6pbAE 
Qt 


(17) 


where 

E  =  AOI  soil  erosion  rate  (m/day) 

Qt  =  total  water  flow  rate  from  AOI  to  all  down-gradient  surface 
waters  (m3/day) 

A  =  AOI  surface  area  (m2) 

pb  =  AOI  soil  dry  bulk  density,  kg/L 

The  factor  1.0E6  in  Equation  17  is  the  conversion  from  kg/L  to  g/m3  or 
mg/L. 

The  equations  above  are  solved  using  the  second  order  accurate  Heun 
method  with  an  adaptive  time  step  of  0.2  days  or  less  to  maintain 
numerical  stability.  The  equation  integration  is  performed  for  daily 
updates  of  both  Crand  TSS. 

Implementation 

Three  AOI  export  treatment  models  have  been  developed  as  stand-alone 
FORTRAN  codes:  degradation  reactor  for  surface  water  and  vadose  zone, 
sedimentation  basin  for  surface  water,  and  sedimentation  basin  and 
reactor  in  tandem  for  surface  water.  The  Tier  2  soil  model  UI  determines 
which  models  to  execute  based  on  user  inputs.  A  processor  within  the  UI 
prepares  input  files  for  each  treatment  model  and  writes  output  files  for 
use  by  down-gradient  models.  This  section  describes  the  processing 
performed  for  implementing  the  aforementioned  AOI  export  treatment 
BMP  models. 

There  are  two  versions  of  the  Tier  2  soil  model,  one  that  uses  average 
annual  hydrology  and  one  that  uses  daily  hydrology.  For  these  reasons, 
mass  and  water  volume  fluxes  can  involve  daily  or  annual  units,  i.e., 
grams/year  or  grams/day  and  cubic  meters/year  or  cubic  meters/day. 

Both  the  reactor  and  sedimentation  basin  models  process  daily  water  flows 
and  mass  fluxes  regardless  of  whether  annual  average  or  daily  hydrology  is 
used.  Thus,  a  soil  model  UI  processor  converts  hydrology  data  into  daily 
units  for  use  in  the  AOI  export  treatment  BMP  models. 
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It  is  possible  that  some  of  the  water  flows  and  mass  fluxes  exiting  the  AOI 
will  not  receive  treatment.  Some  flows  and  fluxes  may  bypass  the 
treatment  reactor  and/or  sedimentation  basin.  Thus,  the  user  is  required 
to  enter  the  fraction  of  total  flow  and  flux  from  AOI  to  surface  water  or 
vadose  zone  that  is  being  treated ,ft.  The  soil  model  UI  processor 
multiplies  the  AOI  total  export  water  flow  rate  and  mass  flux  by  this 
fraction  to  obtain  the  treated  flows  and  fluxes  that  are  input  to  the 
treatment  system(s). 

The  Tier  2  soil  model  generates  five  output  fluxes  to  comprise  AOI  MC  mass 
export  (mass/time)  to  down-gradient  models:  runoff  flux  ( Fr ),  which  is  all 
dissolved;  soil  erosion  flux  (Fe),  which  includes  dissolved  and  particulate 
(soil  adsorbed)  MC;  soil  interflow  flux  (Fj),  which  is  dissolved;  leaching  flux 
(Fi),  which  is  dissolved;  and  the  flux  due  to  erosion  of  solid  phase  MC  (Fes), 
which  is  particulate  since  it  has  not  yet  dissolved.  Soil  erosion  flux  contains 
soil-adsorbed  MCs  and  MCs  dissolved  in  pore  water.  Without  AOI  export 
BMPs,  the  soil  model  UI  processor  accesses  the  five  fluxes  to  create  two 
types  of  output  within  the  water  flux  file  (WFF)  used  by  down-gradient 
models,  one  for  surface  water  and  one  for  the  vadose  zone/groundwater. 
Both  WFF  output  types  contain  MC  mass  fluxes  (grams/year)  and  water 
flows  (cubic  meters/year).  Time  units  of  years  are  used  in  the  WFF  for  both, 
average  annual  and  time-varying,  daily  hydrology.  The  WFF  surface  water 
data  also  contains  two  data  types  for  MC  mass  fluxes,  dissolved  and 
particulate,  whereas  all  of  the  WFF  vadose  MC  mass  fluxes  are  dissolved. 
The  WFF  water  flow  rate  to  vadose  zone,  Qi,  consists  of  net  infiltration  after 
accounting  for  losses  to  ET  and  soil  interflow.  The  WFF  total  water  flow  rate 
to  surface  water  is  the  combination  of  runoff  flow  rate,  Qr,  and  soil  interflow 
flow  rate,  Qi,  which  resurfaces  to  surface  water  flow. 

When  AOI  export  BMPs  are  not  used,  the  soil  model  UI  processor 
combines  Fr  and  Fi  to  form  the  combined  dissolved  flux  within  WFF 
surface  water.  When  AOI  export  BMPs  are  used,  the  soil  model  UI 
processor  excludes  soil  interflow  fluxes  from  those  being  treated  since  it  is 
assumed  that  all  BMPs  will  be  up-gradient  of  return  flow  from  soil 
interflow;  thus,  soil  interflow  is  not  treated  by  AOI  export  BMPs.  Likewise, 
the  total  (treated  and  untreated)  surface  water  flow  rate,  Qt,  excludes 
interflow;  thus,  Qt  =  Qr. 

The  particulate  flux  of  the  WFF  for  surface  water  contains  all  of  the  erosion 
flux,  whereas  a  portion  of  Fe  is  actually  dissolved.  This  misrepresentation 
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presents  no  adverse  consequences  since  the  total  MC  flux  (particulate  and 
dissolved)  is  used  to  compute  reactor  influent  partitioning  and  dissolved 
concentrations,  as  shown  later  below.  Partitioning  of  influent  is  necessary 
since  the  distribution  coefficient  for  surface  water  TSS  can  be  quite  different 
from  that  used  for  soil-pore-water. 

When  an  AOI  export  BMP  includes  a  sedimentation  basin  and/or  a  reactor 
for  surface  water,  a  soil  model  UI  processor  adds  the  fluxes  Fr,  Fe,  and  Fes 
to  form  Fc,  which  in  this  case  is  the  combined  MC  overland  total  flux  for 
surface  water.  The  eroded  solid  phase  MC  overland  flux  Fes  is  included  as  if 
it  is  a  non-solid  flux  although  water  has  not  dissolved  it  at  this  point. 
Assuming  that  water  has  dissolved  it  greatly  reduces  the  complexity  of 
implantation  while  contributing  little  inaccuracy.  Fc  is  simply  the  leaching 
flux  Fi  for  a  vadose  zone  reactor;  of  course,  Fc  is  multiplied  by  ft  to  obtain 
the  treated  flux  that  is  input  to  the  BMP  models. 

The  total  (dissolved  and  particulate)  influent  concentration,  Cn,  (grams/ 
cubic  meter)  required  for  the  sedimentation  basin  model  is  computed  by  the 
soil  model  UI  processor  from  the  ratio  of  AOI  export  flux  and  flow, 

C„=tL  (18, 

where,  as  stated  previously,  Qt  is  the  total  water  flow  rate  from  AOI  to 
surface  water  due  to  overland  runoff.  For  reactors  in  the  vadose  zone,  Qt  is 
simply  the  leaching  flow  rate  Qi  as  stated  previously.  The  reactor  flow  rate 
is  Q  -ft  Qt.  The  untreated  flow  rate  is  (l  -  ft )  QT ,  and  the  untreated  flux  is 

( 1  -  ft) Fc.  The  soil  UI  processor  tracks  untreated  flows  and  fluxes  since 

they  must  be  added  to  the  treated  flows  and  fluxes  in  the  WFF  for  models 
down-gradient  of  the  BMP  models. 

The  dissolved  influent  concentration,  Ci,  for  a  surface  water  reactor  is 
determined  from 


Ct  =  FdCTi 


(19) 


The  fraction  dissolved  Fd  is  computed  from 
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Fd  = 


1  +  10 -6TSSirKdw 


(20) 


where  TSSi,-  is  the  TSS  concentration  (milligrams/liter)  entering  the 
reactor.  Without  tandem  treatment  (i.e.,  no  sedimentation  upstream  of  the 
reactor),  TSSU-  equals  TSSi  as  determined  from  Equation  17.  For  tandem 
treatment,  TSSir  equals  TSS,  which  is  the  computed  TSS  concentration 
exiting  the  sedimentation  basin.  The  latter  distinction  is  handled  within 
the  tandem  model. 


The  soil  model  UI  processor  assembles  and  calculates  the  AOI  water  flow 
rate,  Qt,  and  mass  flux  Fc  for  each  day  and  writes  the  values  to  input  files 
for  the  AOI  export  treatment  models.  Additionally,  the  processor 
calculates  daily  values  of  TSSi  and  writes  these  to  the  same  input  files.  The 
content  of  these  input  files  varies  depending  on  whether  a  reactor, 
sedimentation  basin,  or  tandem  treatment  is  used,  but  there  is  much 
similarity.  Figure  2  shows  an  example  of  the  tandem  treatment  input  file 
produced  by  the  processor. 


Figure  2.  Example  soil  model  UI  processed  input  file  for  AOI  export  tandem 

treatment  model. 


"Test  case  for  generic  sedimentation  basin  and 

reactor  in  tandem  for  range  BMPs" 

"Basin  surf  area,  mA2" 

1000. 

"Basin  mean  depth,  m" 

5. 

"TSS  settling  rate,  m/day" 

2.0 

"Reactor  Length,  m" 

10.0 

"Reactor  Width,  m" 

3.0 

"Reactor  Height,  m" 

1.0 

"Reactor  Porosity" 

0.50 

"Reactor  bulk  density,  kg/L 

1.4 

"fraction  export  treated" 

1. 

"Number  of  MC" 

2 

"MC  name" 

TNT 

"MC  CASID" 

118967 

"TSS-water  Kd" 

1.0 

"Reactor  Kd,  L/kg" 

20.0 

"Reactor  reaction  rate,  1/day" 

10.0 

"Number  of  time  series" 

22280 

Runoff 

AOI 

flow  AOI  flux 

TSS 

year  month  day 

m3/day  g/day 

mg/L 

1950  1  1 

0  0 

16800 

1950  1  2 

0  0 

16800 

1950  1  3 

3224.78  430 

16800 

1950  1  4 

0  0 

16800 

1950  1  5 

0  0 

16800 

1950  1  6 

0  0 

16800 

1950  1  7 

3925.82  430 

16800 

1950  1  8 

0  0 

16800 

1950  1  9 

0  0 

16800 

1950  1  10 

0  0 

16800 

The  soil  model  UI  processor  assembles  and  writes  information  to  the  WFF 
for  use  by  down-gradient  models.  Down-gradient  receiving  water  models 
are  usually  required.  For  example,  there  may  be  a  lake  or  stream  down- 
gradient  of  the  AOI  export  treatment  system  where  MC  concentrations 
must  be  predicted  given  that  not  all  MC  mass  is  removed  by  the  treatment 
system.  The  processor  assembles  and  combines  down-gradient  water  flow 
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rate  and  mass  fluxes  for  each  MC  for  each  day.  Furthermore,  the  processor 
distinguishes  the  mass  fluxes  as  dissolved  and  particulate  as  required  by 
the  WFF  specifications.  The  AOI  export  treatment  models  accommodate 
this  need  as  follows. 

Since  the  reactor  only  removes  dissolved  MCs,  the  portion  of  reactor 
influent  (and  effluent)  that  is  in  particulate  form  (Ftp)  must  be  tracked  so 
that  it  can  be  passed  to  down-gradient  models.  The  reactor  particulate  flux 
is  computed  within  the  reactor  model  from 

F9=(l-Fd)QC„  (21) 

when  a  sedimentation  basin  is  not  upstream  in  tandem.  When  in  tandem 
with  a  sedimentation  basin,  Cn  in  Equation  21  is  replaced  with  Ct,  which  is 
the  outflow  total  concentration  for  a  sedimentation  basin  and  is  available 
within  the  tandem  model. 

Sedimentation  basin  total  mass  flux  out  is  the  product,  QCt.  This  total  flux 
out  is  partitioned  by  the  sedimentation  basin  model  into  particulate  and 
dissolved  fluxes  (FSbP  and  Fsbd,  respectively)  to  maintain  the  two  data  types 
that  are  required  within  the  WFF  for  surface  water.  This  partitioning  is 
computed  as  follows, 


Fsbn  —  FdCtQ 


(22) 


F,u=(l-Ft)CTQ 


(23) 


where  Fp  is  computed  from  Equation  16. 

The  soil  model  UI  processor  uses  the  dissolved  and  particulate  effluent 
fluxes  from  the  AOI  export  treatment  models  to  develop  the  WFF  used  by 
down-gradient  receiving  water  models.  The  treated  output  fluxes  are 
added  to  the  untreated  AOI  export  fluxes  as  well  as  any  soil  interflow 
fluxes,  and  the  combined  results  are  written  to  the  WFF.  The  two  data 
types  for  surface  water  (dissolved  and  particulate)  are  maintained  in  the 
WFF  by  the  processor.  Figure  3  shows  a  schematic  of  the  AOI  export 
treatment  models  and  related  processing. 
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Figure  3.  Schematic  of  AOI  export  treatment  models  and  related  processing. 


Untreated  Flow  and  mass  flux 


WFF  to  Vadose 


As  stated  previously,  three  AOI  export  treatment  models  have  been 
developed  as  FORTRAN  codes:  degradation  reactor  for  surface  water  and 
vadose  zone,  sedimentation  basin  for  surface  water,  and  sedimentation 
basin  and  reactor  in  tandem  for  surface  water.  Each  of  the  three  models 
can  handle  multiple  MCs.  These  three  treatment  models  can  be  applied  in 
seven  possible  combinations: 

1.  Degradation  reactor  for  surface  water 

2.  Degradation  reactor  for  vadose  zone 

3 .  Degradation  reactor  for  surface  water  and  vadose  zone 

4.  Sedimentation  basin  for  surface  water 

5.  Tandem  sedimentation  basin  and  degradation  reactor  for  surface  water 

6.  Tandem  sedimentation  basin  and  degradation  reactor  for  surface  water 
and  vadose  zone  reactor 

7.  Sedimentation  basin  for  surface  water  and  degradation  reactor  for  vadose 
zone 

However,  the  above  combinations  can  be  handled  with  three  possible 
input  options: 

1.  Degradation  reactor  for  surface  water 

2.  Degradation  reactor  for  vadose  zone 

3 .  Sedimentation  basin  for  surface  water 
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If  all  three  options  are  selected,  condition  six  is  invoked.  If  the  first  and 
third  options  are  selected,  condition  five  is  invoked.  If  the  second  and 
third  options  are  selected,  condition  seven  is  invoked,  and  so  on. 
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3  Inputs  and  outputs 

Input  and  output  for  the  TREECS™  BMP  modules  are  presented  in  this 
chapter.  This  information  is  organized  in  the  same  fashion  as  done  in 
Chapter  2,  source  removal  and  AOI  export  treatment  BMPs. 

Source  removal  BMP  inputs 

The  input  information  associated  with  each  of  the  five  source  removal 
BMPs  is  presented  in  this  section. 

Soil  removal  and  landscape  burning 

The  Tier  2  soil  model  UI  has  a  new  tab  labeled  Source  Removal  BMPs.  The 
inputs  for  all  five  source  removal  BMPs  are  entered  via  this  tab.  There  are 
four  sub-tabs  on  the  Source  Removal  BMPs  screen.  These  are  labeled  Soil 
Removal/ Burning,  Phytotransformation,  Phytoextraction,  and  Selective 
MC  Removal  of  Chunks.  Figure  4  shows  a  preliminary  version  of  the  input 
screen  with  selection  tabs.  This  version  has  phytotransformation  as  a 
separate  tab  rather  than  as  a  sub-tab  under  Source  Removal  BMPs.  For 
the  final  software  version,  the  phytotransformation  BMP  will  be  moved  to 
a  sub-tab. 

Soil  removal  and  burning  of  the  landscape  are  combined  on  the  same 
screen  since  the  input  requirements  are  conducive  to  combining  them.  For 
the  final  software  version,  the  user  will  be  required  to  check  whether  to  set 
Rns  =  Rs  or  set  Rns  to  zero  for  soil  removal.  This  feature  was  not  included  in 
the  preliminary  input  screen  shown  in  Figure  4.  Setting  Rns  =  Rs  means  the 
soil  is  permanently  removed.  Setting  Rns  to  zero  means  that  solid  MC 
particles  are  removed,  and  the  soil  is  put  back  in  place;  thus,  non-solid 
phase  MCs  are  not  removed. 

Each  MC  included  in  the  application  is  displayed  along  with  its  CASRN.  As 
shown  in  Figure  4,  the  user  selects  the  MC  to  which  burning  applies.  In  the 
lower  panel  of  the  screen,  the  user  enters  the  soil  removal  rates  in  metric 
tons  per  year  and  the  burning  rate  in  acres  per  year  for  each  year  that 
there  are  changes  in  these  two  inputs.  At  least  two  time  points  in  years 
must  be  entered,  and  inputs  are  stepped  for  each  input  year,  i.e.,  they  are 
constant  until  the  next  update  year.  The  UI  displays  the  computed  values 
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for  Rs  and  Rns  for  soil  removal  and  burning  and  checks  to  make  sure  the 
calculated  values  and  their  combination  for  soil  removal  and  burning  do 
not  exceed  1.0  for  any  input  year. 

Figure  4.  Preliminary  input  screen  for  soil  removal  and  landscape  burning  BMPs. 

Site  Characteristics  |  Soil  Properties  Hydrology  Fate/Transport  Parameters  [  Chemical-Specific  Properties  \  Miscellaneous  BMP  Info  AOI  Phytotransformation  BMP  Source  Removal  BMPs 

Soil  Removal/Buming  [  Phytoextraction  j  Selective  MC  Removal  of  Chunks 


Select  MC  to  which  burning  applies  (if  applicable): 


►* 
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Phytotransformation 

The  input  screen  for  phytotransformation  displays  each  MC  included  in  the 
application  along  with  its  CASRN.  The  user  selects  each  MC  that  is  treated 
with  phytotransformation  and  MC-specific  values  for  the  plant  growth  rate 
G,  bioconcentration  ration  BCR,  and  fraction  of  MC  mass  taken  up  by  the 
plant  that  is  transformed, /p.  Each  of  these  three  inputs  is  constant  over 
time  for  each  MC.  Although  not  shown  in  the  phytotransformation 
preliminary  input  screen  of  Figure  5,  an  input  table  is  provided  for  entering 
the  time-varying  fraction  of  AOI  area  that  is  treated  with  phytotransforma - 
tion,yp.  This  table  will  be  similar  to  the  one  described  above  for  soil  removal 
and  burning.  After  entering  the  time  pairs  for/p ,  the  user  can  view  the 
computed  MC  removal  rate  and  half  life  due  to  phytotransformation  for 
each  time  point.  The  UI  checks  to  ensure  that  none  of  th e/p  inputs  for  any 
MC  is  greater  than  1.0  for  all  time  pairs.  Additionally,  the  fraction/pis  not 
allowed  to  exceed  1.0. 
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Figure  5.  Preliminary  input  screen  for  phytotransformation  BMP. 


Model  Options  _ 

Site  Characteristics  ]  Soil  Properties  ]  Hydrology  ]  Fate/Transport  Parameters  J  Chemical-Specific  Properties  [  Miscellaneous  BMP  Info  |  AO  I  Phytotransformation  BMP  Source  Removal  BMPs  |  ADI  Export  Treatment  BMPs 
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Phytoextraction 

As  for  the  input  screen  described  above,  the  input  screen  for 
phytoextraction  displays  each  MC  included  in  the  application  along  with 
its  CASRN.  The  user  selects  each  MC  that  is  treated  with  phytoextraction 
and  enters  the  plant  production  rate  and  bioconcentration  ratio  for  that 
MC.  An  input  table  of  time  pairs  for  plant  harvesting  fractional/  is 
provided  for  each  selected  MC,  where  the  user  enters  the  year  and  the 
fraction.  As  explained  in  Chapter  2,  the  plant  harvesting  fraction  is 
actually  the  fraction  of  AOI  area  that  is  treated  with  phytoextraction.  As  a 
result,  the  UI  checks  to  ensure  that  none  of  the/#  inputs  for  any  MC  is 
greater  than  1.0  for  all  time  pairs.  Figure  6  shows  the  preliminary  input 
screen  for  this  BMP.  The  user  can  also  view  the  calculated  phytoextraction 
removal  rate  Rns  for  each  input  year. 


Figure  6.  Preliminary  input  screens  for  phytoextraction  BMP. 


Selective  MC  removal 

Again,  this  input  screen  displays  each  MC  included  in  the  application 
along  with  its  CASRN,  and  the  user  selects  each  MC  that  is  selectively 
removed  via  removal  of  mass  chunks.  For  each  selected  MC,  there  is  a  time 
pair  input  table  for  the  user  to  enter  the  year  and  the  removal  rate  SR 
(grams/year)  as  shown  in  Figure  7. 
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Figure  7.  Preliminary  input  screens  for  selective  MC  removal  BMP. 


Site  Characteristics  j  Soil  Properties  [  Hydrology  |  Fate/Transport  Parameter  |  Chemtcal-Speaftc  Properties  j  Miscellaneous  BMP  Into  |  Hjl  PhytotransfonTiation  BMp]  Source  Removal  BMPs  j 

Soil  Removal/Buming  |  Phytoextraction  !  Selective  MC  Removal  of  Chunks  | _ 


Note:  Removal  of  duds  can  be  handled  within  the  ‘Site  Conditions/Operational  Inputs'  screen  by  setting  percent  sympathetic  duds  to  0  within  the  'Munitions  usage  information'  panel. 
Select  MC  to  which  selective  removal  applies  and  provide  rates: 
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AOI  export  treatment  BMP  inputs 

The  input  information  associated  with  the  AOI  export  treatment  BMPs  is 
presented  in  this  section. 

Degradation  reactor 

Input  for  the  reactor  model  include  the  following  reactor  parameters: 
length  Lr  (meters),  width  Wr  (meters),  height  Hr  (meters),  porosity  dry 
bulk  density  pbr  (kilograms/liter),  and  fraction  of  AOI  surface  water  export 
that  is  treated,  ft.  Three  MC-specific  inputs  are  required,  which  include: 
reactor  sorption  distribution  coefficient  Kdr  (liters /kilogram)  for 
partitioning  of  MC  between  pore  water  and  reactor  solid  medium;  reactor 
reaction  rate  Ar  (per  day);  and  sorption  distribution  coefficient  Kdw 
(liters/kilogram)  for  partitioning  of  MCs  between  water  and  TSS.  Figure  8 
shows  the  preliminary  input  screen.  The  input  for  Kdw  is  not  needed  or 
used  for  a  vadose  zone  reactor,  but  this  input  field  is  retained  since  the 
same  UI  and  model  are  used  for  surface  water  and  vadose  zone  reactors. 
The  values  for  the  MC  sorption  partitioning  coefficients  Kdr  and  Kdw  can  be 
quite  different.  Likewise,  values  of  Kdw  can  be  different  from  partitioning 
coefficients  for  soil  and  sediment  pore  water. 
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Figure  8.  Preliminary  input  screen  for  degradation  reactor  BMP. 
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Sedimentation  basin 

Input  for  the  sedimentation  basin  model  include  the  following  basin 
parameters:  surface  area  Ab  (square  meters),  mean  depth  Hb  (meters),  TSS 
settling  rate  vs  (meters/day),  and  fraction  of  AOI  surface  water  export  that 
is  treated  ft.  The  basin  volume  Vb  used  in  the  model  is  the  product  AbHb. 
One  MC-specific  input  is  required,  which  is  the  sorption  distribution 
coefficient  Kdw  (liters /kilogram)  for  partitioning  of  MCs  between  water 
and  TSS.  Figure  9  shows  the  preliminary  input  screen. 


Figure  9.  Preliminary  input  screen  for  sedimentation  basin  BMP. 


Vadose  zone  reactor  and  tandem  treatment 

The  inputs  for  the  vadose  reactor  are  identical  to  the  reactor  inputs 
presented  above  except  for  the  labeling,  i.e.,  surface  water  (SW)  versus 
vadose  zone  (VZ).  The  inputs  for  the  tandem  model  are  merely  a 
combination  of  inputs  for  a  reactor  and  sedimentation  basin  except  that 
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the  input  for  fraction  treated  by  the  reactor  is  eliminated  since  it  is  a 
redundant  input. 

Outputs 

Source  removal 

Source  removal  BMPs  affect  AOI  soil  MC  mass  and  concentration  and  MC 
mass  transported  from  the  AOI;  thus,  there  are  no  direct  outputs  associated 
with  source  removal  BMPs.  However,  there  are  the  inputs  calculated  by  the 
UI  for  input  to  the  soil  model  that  were  described  previously  in  this  chapter. 
The  effects  of  source  removal  BMPs  are  manifested  by  changes  in  AOI  soil 
MC  mass  and  concentration  and  mass  flux  exiting  the  AOI. 

AOI  export  treatment 

Unlike  the  source  removal  BMPs,  each  of  the  AOI  export  treatment  BMPs 
has  viewable  output.  This  output  is  produced  by  the  three  AOI  export 
treatment  models  (reactor,  sedimentation  basin,  and  tandem  sedimentation 
basin,  followed  by  reactor),  and  it  is  used  to  process  input  via  WFFs  for  any 
down-gradient  receiving  water  model.  Users  can  view  AOI  export  treatment 
output  by  going  to  the  View  Results  tab  and  selecting  AOI  Export  treatment 
BMP  Output  under  the  Soil  Viewer  of  TREECS™.  Figure  10  shows  an 
example  output  for  tandem  treatment.  The  output  includes  an  echo  print  of 
input  parameters  and  tabular  daily  values  for  each  MC  of:  entering  MC 
mass  flux  (Auxin),  entering  MC  total  concentration,  MC  total  concentration 
exiting  the  sedimentation  basin,  MC  total  concentration  exiting  the  tandem 
system,  MC  mass  Aux  exiting  the  tandem  system,  particulate  (part)  mass 
Aux  delivered  to  down-gradient  (DG)  models,  dissolved  (diss)  mass  Aux 
delivered  to  DG  models,  TSS  concentration  in  the  sedimentation  basin,  and 
sedimentation  basin  time  step.  Output  for  the  reactor  and  sedimentation 
basin  are  similar  but  involve  less  information  than  that  produced  by  the 
tandem  model.  Output  of  all  three  models  includes  down-gradient  dissolved 
and  particulate  Auxes. 
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Figure  10.  Example  output  for  tandem  AOI  export  treatment  BMP. 
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4  Summary 

This  report  provides  formulations  used  for  the  source  removal  and  AOI 
export  treatment  BMPs  and  guidelines  for  their  final  implementation  within 
TREECS™.  Work  was  initiated  but  not  completed  on  the  coding  to 
implement  these  BMPs  within  TREECS™.  Following  the  initial 
implementation  within  TREECS™,  some  of  the  formulations  were  modified, 
and  improvements  were  made  in  the  plans  for  interfacing  these  BMP 
modules  within  TREECS™.  Thus,  this  report  serves  to  document  the  final 
formulations  and  provide  information  for  final  implementation  within 
TREECS™.  This  report  also  documents  the  BMP  module  input 
requirements  and  output. 
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